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Let S be a regular set of R'' and X : S — >■ M be Gaussian field with regular paths. In order to give 
bound to the tail of the distribution of the maximum, we use the record method of Mercadier. We 
present some new form in dimension 2 and extend it to dimension 3 using the result of the expectation 
of the absolute value of quadratic forms by Li and Wei. Comparison with other methods is conducted. 
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1 Introduction 



The problem of computing the tail of the maximum has a lot of applications in spatial statistics, 
image processing, oceanography, genetics etc see for example Cressie and Wikle [7]. It is exactly 
^— ( solved only for about ten processes with parameter of dimension 1, see Azais and Wschebor [5] p4 

^ for a complete list. In the other cases, one has to use some approximations. Several methods have 

been used, in particular 

^ • The tube method. Sun [IT]- 

• Double sum method, Piterbarg [14) . 

• Euler characteristic method see, for example, Adler and Taylor [I]. 

• Rice or direct method, Azai's and Delmas [3], Azais and Wschebor [5]. 

With respect to these methods, the record method which is the main subject of this paper and which 
is detailed in Section [2] has the advantage of simplicity and also the advantage of giving a bound 
which is non asymptotic: it is true for every level and not for large u only. 

It has been introduced for one-parameter random processes by Rychlik and extended to two- 
parameter random fields by Mercadier |13j to study the tail of the maximum of smooth Gaussian 
random fields on rather regular sets. 

It has two version, one is an exact implicit formula : Theorem 2 in |13| that is interesting for numerical 
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purpose and that will not be considered here; the other form is a bound for the tail, see inequality 
([T| hereunder. 

This bound has the advantage of its simplicity. In particular it avoids the computation of the 
expectation of the absolute value of the Ifessian determinant as in the direct method of but it 
works only dimension 2. 

For practical applications, the dimensions 2 and 3 (for the parameter set) are the most relevant so 
there is a need of an extension to dimension 3 and this is done in Section 3 using results on quadratic 
forms by Li and Wei 

The bound also has the drawback of demanding a parameterization of the boundary. For example, if 
we consider the version of Azai's and Wschebor (:5 , Theorem 9.5 ) of the result of Mercadier, under 
some mild conditions on the set S' C and on the Gaussian process X, we have 

P{Ms > u} <P{Y{0) >u}+ f E(jy'(/)| I Y{1) = u)py(i){u) dl 

Jo 

+ f E(|X('i(t)-X^(t)+l I X{t)=u, = 0)pxw,x;(t)(«,0)di, (1) 

Js 

where 

• Ms is the maximum of X{t) on the set S. 

• Y{1) = X{p{l)) with p : [0, L] — >■ dS is a parameterization of the boundary dS by its length. 
^„ _ d'X 



dxidxj ' 

• pz{x): the value of the density function of random vector Z at point x. 

• x"*" = sup(a;, 0), = sup(— x, 0). 

The proof is based on considering the point with minimal ordinate (second coordinate) on the level 
curve. As we will see, this point can be considered as a "record point". 

So the second direction of generalizations is to propose nicer and stronger forms of the inequality ([T]). 
This is done in Section [2] The result on quadratic form is presented in Section |4] and some numerical 
experiment is presented in Section [5] 

Notation 

o 

• S is some rather regular set included in or R'^. dS is its boundary; S is its interior. 

• Ms ~ maxX(s) where X{s) is some rather regular process. 

s£S 

• CTi is the surface measure of dimension i. It can be defined as a Hausdorff measure. 

• X' , X" are the first and second derivatives of the process X{t). In particular if a is some 
direction then X'^ is the derivative along the direction a. 

• M ^ means that the square matrix A/ is semi-definite negative. 

• S'^'^ is the tube around S, i.e 

= {s e R^ : dist(s,S') < e}. 

• dn is the Hausdorff distance between sets, defined by 

dniS, T) = inf{e : S C T+^ T C S+'j. 

• ip{x) and <^{x) are the density and distribution function of a standard normal variable. 
^{x) = l-^{x). 
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2 The record method in dimension 2 revisited 



We will work essentially under the following assumption: 
Assumption 1: {X{t), t £ NS C R^} is a Gaussian stationary field, defined in a neighborhood NS 
of S with paths and such that there exists some direction, that will be assumed (without loss of 
generality) to be the direction of the first coordinate, in which the second derivative Xi'i{t) exists. 
We assume moreover the following normalizing conditions that can always be obtained by a scaling 

E{X{t)) = 0, Var(X(t)) = 1, VarX'(t) = h. 

Finally we assume that Var(X"i(t)) > 1 which is true as soon as the spectral measure of the process 
restricted to the first axis is not concentrated on two opposite atoms. 

In some cases we will assume in addition 
Assumption 2: X{t) is isotropic, i.e Cov{X{s), X{t)) = p{\\t — s\\^), with paths and S is a convex 
polygon. 

Under Assumption 1 and 2 plus some light additional hypotheses, the Euler Characteritic (EC) 
method [T] gives 

P{Ms >u} = Psiu) + Rest, 

with 

2\J2-K 27r 

where the rest is super exponentially small. 
The direct method gives [4] 

P{Ms > «} < Vm(u) = $(w) + ^li^ / \cip(xlc) + x^[xlc)\ <p{x)dx 

2\/2tv Ju 



2n 



'(0))^/^exp(-:c^(24p"(0)-2)-^) 
v/24p"(0) - 2 



ifi{x)dx, (2) 



where c = ^Var(X('i) - 1 = ^12p"(0) - 1. 
The record method gives l.'i 

P{Ms >u}< $(tt) + "^^^ifiiu) + IMu/c) + u^u/c)] ^{u). 

V2tt Stt 

A careful examination of these equations shows that the main terms are almost the same except that 
in the record method the coefficient of ai{dS) is twice too large. When 5* is a rectangle [0, Ti] x [0, T2], 
it is easy to prove that this coefficient 2 can be removed, see for example Exercise 9.2 in [5]. 

The goal of this section is to extend the result above to more general sets and to fields satisfying 
Assumption 1 only. The main result of this section is the following 

Theorem 1. Let X satisfy the Assumption 1 and suppose that S is the Hausdorff limit of connected 
polygons S„. Then, 

P{Ms >u}< $(.) + ^^^■^^^nO.idS^Mu) ^ ^^^^^^^ ^ ^^^^^ 

2v27r 27r 

where c = VVar(X('J - 1. 

Remark: the choice of the direction of ordinates is arbitrary and is a consequence of the arbitrary 
choice of the the second derivative When the process X{t) admits derivative in all direction, 

the choice that gives the sharpest bound consists in chosing as first axis, the direction a such that 
Var(X^'Q,) is minimum. 
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Unfortunately the proof it is based on an exotic topological property of the set S that will be 
called "emptyable". 

Definition 1. The compact set S is emptyable if there exists a point O £ S which has minimal 
ordinate, and such that for every s £ S there exists a continuous path inside S from O to s with non 
decreasing ordinate. 

In other word, suppose that 5* is filled with water and that gravity is in the usual direction; S is 
emptyable if after making a small hole at O, all the water will empty out, see Figure [l] 



Figure 1: Example of non-emptyable set. The non-emptyable part is displayed in black. 



Proof. Stepf : Suppose for the moment that X has C°° paths and that S is an emptyable polygon. 
Considering the event {AIs > u}, we have 

P{Ms >u}^ P{X{0) >u}+ P{X{0) < u, Ms > u}. (4) 

It is clear that if X{0) < u and Ms > u, because S is connected, the level curve 

Cu = {t£S: X{t) = u} 

is not empty, and there is at least one point T on Cu with minimal ordinate. There are two possibilities: 

• T is in the interior of S. In that case, suppose that there exists a point s £ S with smaller 
ordinate than T (s2 < T2), such that X{s) > u. Then, due to the emptyable property, on the 
continuous path from O to s there would exist one point s' with smaller ordinate than T, and 
with X{s') — u. This is in contradiction with the definition of T. So we have proved that for 
every s £ S , S2 < T2 we have X{s) < u. It is in the sense that T can be considered as a record 
point. It implies that 

{X[{T) = 0, X'^iT) > 0, XUT) < 0}. 
The probability that there exists such a point is clearly bounded, by the Markov inequality, by 

E (card{t G S : X{t) = u, X[{t) = 0, X'iit) > 0, X[\{t) < 0}) . 

Applying the Rice formula to the field Z — {X,Xi) from to R^, we get that 

o 

P{3t £ S : X(t) = u, t has minimal ordinate on Cu} 
< y^E(|det(Z'(t))|Ix,(,)>oI^»(,)<o I = (u,0)) xpz(t)(u,0) 

= a2iS)^ E {X'^,-{t)X'2+{t) \X{t) = u,X[{t) = 0) 

V An 

= a2iS)^ E {X^+{t)) E (Xrr W I X{t) = u,X[{t) = O) 

V An 

= a2{S)^[c^(u/c)+u^u/c)]. (5) 
An 
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Note that the vahdity of the Rice formula holds true because the paths are of class C°° and 
that X{t) and X'i{t) are independent. The computations above use some extra independences 
that are a consequence of the normalization of the process. The main point is that, under the 
conditioning 

det{Z'{t))=X'Mt)XUt). 

T is on the boundary of S that is the union of the edges {Fi, . . . ,Fn)- It is with probability 1 
not located on a vertex. Suppose that, without loss of generality, it belongs to Fi. Using the 
reasoning we have done in the preceding case, because of the emptyable property, it is easy to 
see that 

{X(T) = u, XUT)>0, X'^{T)<0}, 

where a is the upward direction on Fi and /3 is the inward horizontal direction. Then, apply 
the Markov inequality and Rice formula in the edge Fi, 

P{3 1 £ Fi : X{t) = u, t has minimal ordinate on Cu} 

< P{3teFi: X{t) = u, XUt) > 0, X'p{t) < 0} 

< E (card{f G Fi : X{t) = u, X'^{t) > 0, X'p{t) < 0}) 

= y" E (|X(i))!Ixi(t)>oIx^(t)<o 1 = W) X Pxit){u) dt 
L(FiV(n)E(x;,+ (t)I^,(,)<o). 



Denote by 9\ the angle (q,/3). X'p can be expressed as 

cos 6*1 X'^ + sin 9i Y, 
with y is a standard normal variable that is independent with X'^. Then 

E(X;+(t)I;,,(,)<0) 

= ^{X'^ \cos0iX'^+siue-iY<o) 

1 — COS 6\ 

Summing up, the term corresponding to the boundary of S is at most equal to 

^ {1- COS ei)cji{Fi) ip{u)(Ji{dS) 
'^("^^ 2jr. = ' 



since cosOi is just the length of the oriented projection of the boundary of S on the 



X -axis, so it is zero. 



Hence, summing up ([sjl, (|6| and substituting into Q, we obtain the desired upper-bound in our 
particular case. 

Step 2: Suppose now that S is a general connected polygon such that the vertex O with minimal 
ordinate is unique. We define Si as the maximal emptyable subset of S that contains O. It is easy to 
prove that Si is still a polygon with some horizontal edges and that S\Si consists of several polygons 
with horizontal edges, say S2 , . . . , Sj ^ , see Figure [2] 

So we write 

P{Ms >u}< P{X(0) > It} + P{A/si > u, X{0) <u} + ^P{Msi < u, Mg, > u}. (7) 
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Figure 2: Example on construction of Si. 



Suppose for the moment that all the SI, i ~ 1, . . . ,n are emptyable. Then, to give bounds to the 
event 

{Msi < u, M5, > u}, 

we can apply the reasoning of the preceding proof but inverting the direction: in S2, we search points 
on the level curve with maximum ordinate. Let E be the common edge of Si and S^- Clearly, 
when {Msi < u, Mgi > u}, the level curve is non empty and by the same arguments as in Theorem 
[1] there exists t £ S2 satisfying whether (except events with zero probability) 

• t is in the interior of 5*2 and 



{X{t) = u, X[{t) = 0, X'2{t) < 0, x'Mt) < 0}. 
From Markov inequality and Rice formula, this probability is at most equal to 

^(^>f^^ [c^iu/c) + uHu/c)]. (8) 

t lies on some edges of 5*2. Note that t can not belong to E. Then, as in Theorem [l] we consider 
the event t is on each edge and sum up the bounds to obtain 

P ({3i € 882 '■ X{t) = It, thas maximal second ordinate on the level curve} n {Msi < m}) 



2V27r 

(9) 

From in and Q we have 



Summing up all the bounds as in ( |10[ ), considering the upper bound for P{X{0) < u, Ms^ > u} as 
in Theorem [1] and substituting into |7|, we get the result. 

In the general case, when some 5*2 is not emptyable, we can decompose SI as we did for S and by 
induction. Since the number of vertices is decreasing, we get the result. 

Step 3: Passing to the limit. The extension to process with non C°° paths is direct by an 
approximation argument. Let X^it) be the Gaussian field obtained by convolution of X{t) with a 
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size e convolution kernel (for example a Gaussian density with variance We can apply the 

preceding bound to the process 



y^Var(X,(t)) 



where Eg = Ya.r{X ^(t)). Since Var(Xe(t)) — >■ 1 and — > /2 , ma-xt^s Xc{t) — >■ Ms and we are done. 
The passage to the limit for Sn tending to S is direct. □ 

Some examples 

• If 5* is compact convex with non-empty interior then it is easy to construct a sequence of polygons 
Sn converging to S and such that liminf„ (7i(9S'n) = ui{dS), giving 

P{Ms >u}< Pr{u) = + ^l^^{u) + Mu/c) + u'i>{u/c)] ^(u). (11) 

2v27r 27r 

• More generaly, if 5* is compact and has a boundary that is piecewise-C^ except for a finite 



number of points and the closure of the interior of S equals to S, we get (111 by the same tools 



• Let us now get rid of the condition S — S but still assuming the piecewise-C^ condition. Define 
the "outer Minkowski content" of a closed subset 5 C as (see 

OMC(S) = lim^^^^^l\^, 

whenever the limit exists (for more treatment in this subject, see [2]). This definition of the 
perimeter differs from the quantity ai (dS) . A simple counter-example is a set corresponding to 
the preceding example with some "whisker" added. Using approximation by polygons, we get 

P{Ms >u}< Pr{u) = $(i*) + QMC(S) ^^^^ ^ a2{S) ^ ^t$(^^/c)] ^(u). (12) 

2v27r 2-K 

• The next generalization concerns compact r-convex sets with a positive r in the sense of [8]. 
These sets satisfy 

5= Pi m.^\B{x,r). 

S(a;,r)nS = 

This condition is slightly more general than the condition of having positive reach in the sense of 
Federer [9]. Suppose in addition that S satisfies the interior local connectivity property: there 
exists ao > such that for all < a < ao and for all x £ S, int {B{x, a) n S) is a non-empty 
connected set. Then we can construct a sequence of approximating polygons in the following 
way. 

Let Xi,X2, . . . , Xn be a random sample drawn from a uniform distribution on 5* and Sn be the 
r-convex hull of this sample, i.e 

Sn= n M'\B(a;,r), 

B{x,r)n{Xi,X2,...,Xn} = 9 

which can be approximated by polygons with an arbitrary error. By Theorem 6 of Cuevas et 
al [S], S'n is a fully consistent estimator of S, it means that dniSn, S) and dnidSn, dS) tend to 
as n tends to infinity. This implies a2(Sn) — > o-2(S') and OMC(Sn) — ^ OMC(S). Hence, we 



obtain ( 12 1 
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A complicated case: a "Swiss cheese". Here, we consider an unit square and inside it, we 

oo 

remove a sequence of disjoint disks of radius n such that ti"^^''? < 1 to obtain the set S. When 

i=l 

oo 

y^ri < oo the bound ( 3 1 makes sense directly. But examples can be constructed from the 



Sierpinski carpet (see Figure 



3 1 such that = oo : divide the square into 9 subsquares 

of the same size and instead of removing the central square, remove the disk inscribed in this 
square and do the same procedure for the remaining 8 subsquares, ad infinitum. 



Figure 3: Sierpinski carpet (source: Wikipedia). 



In our case, 

oo 



1 V — - - 

4 32i - 4 

i = l 1=1 



This proves that the obtained set S has positive Lebesgue measure and is not fractal. We have 
on the other hand 

>^ ri = - >^ — — = oo. 

i — l i — l 

Let Sn be the set obtained after removing the n-th disk. Since S C Sn, 

, ^ n n 

P{Ms >-!*}< P{Ms„ >u}< $(?i) + ^^(4+27r^ri) + (l-7r^r-?) [cip{u/c) + it$(u/c)] ^{u)/{2tv) 
Hence, 



P{Ms >u}< $(u)+min 



Remarks: 



2v^ 



(4 + 27r^ri) + (1 - vr^r,^) [af{u/c) + u${u/c)] ifi(u)/{2n) 



1. In comparison with other results, all the examples considered here are new. Firstly the condi- 
tions on the process are minimal and weaker than the ones of the other methods. Secondly the 
considered sets are not covered by any other methods. Even for the first example, because we 
do not assume that the number of irregular points is finite, which is needed, for example, for 
the convex set to be a stratified manifold as in [T]. 

2. Theorem [l] can be extended directly to non connected sets using sub-additivity 

P{Ms,uS2 >u}< P{Msi >u} + P{Ms2 > u}. 
This implies that the coefficient of "I>(it) in |3| must be the number of components. 
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Is the bound sharp? 

• Under Assumption 2, Adler and Taylor QJ show that 

liminf - 2u'^log\P{Ms > n} - Pe{u)\ > 1 + 

u—>-\-oo 

From 

< Pfi(u) - Pe{u) = ^^Vi^) [Mu/c) - u^{u/c)] 
and the elementary inequality for a; > 0, 



it is easy to see that 

u— >- + oo 



liminf - 2ulog(PR{u) - Pb(m)) > 1 + 1/c 

So the upper bound P_r(u) is as sharp as P_e(m) . 

• Let S be a compact and simply connected domain in having a C^-piecewise boundary. Assume 
that all the discontinuity point are convex, in the sense that if we parametrize the boundary in 
the direction of positive rotation, then at each discontinuity point, the angle of the tangent has 
a positive discontinuity. Then, it is easy to see that the quantity 

, . dist(s-t,Ct) 

K(b) = sup sup 77^ j— — - 

is finite, where dist is the Euclidean distance and Ct is the cone generated by the set of directions 



A G : |1A|1 = 1, 3s„ G S such that Sn t and |— 1- A 

1 1 Sn 1 1 

In order to apply the Theorem 8.12 in 15', besides the Assumption 1, we make some additional 
assumptions on the field X such that it satisfies the conditions (A1)-(A5) page 185 in [5]. 
Assume that 

— X has paths. 

— The covariance function r{t) satisfies \r{t)\ 7^ 1 for all t ^ 0. 

— For all s ^t, the distribution of {X{s), X{t), X' {s), X' (t)) does not degenerate. 
With these hypotheses, we can see that 

— The conditions (A1)-(A3) are easily verified. 

— The condition (A4) which states that the maximum is attained at a single point, can 
be deduced from Proposition 6.11 in |5] since for s 7^ t, {X{s), X{t), X' (s), X' (t)) has a 
nondegenerate distribution. 

— The condition (A5) which states that almost surely there is no point t £ S such that 
X'{t) — and det{X" (t)) — 0, can be deduced from Proposition 6.5 in [5] applied to the 
process X'{t). 

Since all the required conditions are met, by Theorem 8.12 in [51, we have 

liminf - 2x'^ log \Pm{x) - P{Ms > x}] > 1 + inf „ ^ „ > 1, (13) 

where 

a?= sup Var(X(s)|X(t),X'(t)) 



ses\{t} (l-^(s,i)) 



2 
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and 

dist(|:r(s,t),a) 

Kt = sup , " I ,s -■ 



Note that the condition ^(S) is finite implies that /t(t) is also finite for every t £ S. ( 13 1 is true 
also for Pr, since as x — >■ +00, Pr{x) is smaller than Pm{x) (see Section [s] for the easy proof). 
As a consequence P_r is super exponentially sharp. 

• Suppose that 5 is a circle in R^. Then {X{t) : t G 5} can be viewed as a periodic process on 
the line. In that case, it is easy to show, see for example Exercise 4.2 in [5], that as it 00 

P{Ms > u) = ^^(u) + 0(^(n(l + S)) = ^^^^(u) + OMuil + S)) 



for some 5 > 0; while Theorem [T] gives with a standard approximation of the circle by polygons 

P(Ms >u)< Pn(u) = $(ii) + 2M2M^(^), 

2\/27r 



which is too large. This shows that the bound Pr is not always super exponentially sharp. 



3 The record method in dimension 3 

For example, with the direct method, some difficulties arise in dimension 3 because we need to 
compute 

E| det(X"(t)|, 

under some conditional law. This can be conducted only in the isotropic case using random matrices 
theory, see [4| and even in this case the result is complicated. In dimension 2, the record method is 
a trick that permits to spare a dimension in the size of the determinant we have to consider because 
the conditioning implies a factorization. For example in equation ([5| we have used the fact that 

Aet{Z'(t)) = X'i^(t)X'^{i), 

under the condition. In this section we will use the same kind of trick to pass from a 3,3 matrix 
to a 2,2 matrix and then a 2,2 determinant is just a quadratic form so we can use, to compute the 
expectation of its absolute value, the Fourier method of Berry and Dennis |6| or Li and Wei [llj . 
This computation is detailed in Section 4 and is one of the main contributions of the paper. 
Before stating the main theorem of this section, we recall the following lemma (see Chapter 5 of 
Prasolov and Sharygin [15) ) 

Lemma 1. Let Oxyz be a trihedral. Denote by a, b and c the plane angles xOy, yOz and zOx, 
respectively. Denote by A, B and C the angles between two faces containing the line Oz, Ox and Oy, 
respectively. Then, 

a. sin a : sin A — sin 6 : sin_B = sine : sinC . 

b. cos a = cos bcos c + sin 6 sin ccos j4. 

Our main result is the following 

Theorem 2. Let S be a compact and convex subset o/R^ with non-empty interior and let X satisfy 
Assumption 1. Suppose, in addition that X is isotropic with respect to the first and second coordinate, 
i.e 

Cov{X{t-i,t2,t3); X{si,S2,t3)) = p{{ti — Si)^ + {t2 — 82)^) with p of class . 
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Then, for every real u, 



V{M >u]< 1 - 



/27r 



47r 



v/l2p"(0) - 



+ 



Cr3(5')(/3(u) 



(27r)3/2 

where A is f/ie caliper diameter. 



1 + 



(8p"(0))^/^exp(-^.^(24p"( 0) - 2)'^) 
V24p"(0) - 2 



Vl2p" 


(o)-i 


2)")] 





Vl2p"(0)-1 



Proof. By the same limit argument as in Theorem [T] we can assume that X[t) has C°° paths and 
that 5 is a convex polyhedron. Let O be the vertex of 5* that has minimal third coordinate, we can 
assume also that this vertex is unique. It is clear that if X{0) < u and Ms > u then the level set 

C{u) = {te S: X{t) = u} 

is non empty and there exists at least one point T having minimal third coordinate on this set. Then, 

P{Ms >u} = P{X{0) >u} + P{X{0) <u, Ms> u} 

< P{X{0) > u} + P{3T e S : X(T) = u, Thas minimal third coordinate onC„}. 

(14) 

Now, we consider three possibilities: 

• Firstly, if T is in the interior of S, then by the same arguments as in Theorem [T] for all the point 
s £ S with the third coordinate smaller than the one of T, X{s) < X{T); it means that, at T, X{t) 
has a local maximum with respect to the first and second coordinates and is non-decreasing with 
respect to the third coordinate. Therefore, setting 



A{t) = 



x'Mt) xut) 
x'Mt) xUt) 



we have 

{X{T) = u, X[iT) = 0, X^iT) = 0, A{T) ^ 0, X^(T) > 0}. 
Then, apply the Rice formula to the field Z = {X, X[, X2) and the Markov inequality, 

o 

P{3r G S : X{T) = u, Thas minimal third coordinate onC„} 

< P{3t G 5 : X{t) = It, X[{t) = 0, X2{t) = 0, Xs{t) > 0, A{t) < 0} 

< E(card{f G 5" : X{t) = u, X[{t) = 0, X'2{t) = 0, X^{t) > 0, A{t) ^ 0}) 
= E(card{f G 5" : Z{t) = {u, 0, 0), X'^it) > 0, A{t) < 0}) 

E(\det{Z'{t))\Ixut)<o^A(t)^o \ Z{t) = (m,0,0)) x pz(t)(u,Q,G) dt. 



Under the condition Z{t) = (ii,0,0), it is clear that det(Z'(t)) = X'2,{t) Aet{A(t)). So, we obtain the 
bound 

oz{S) ^ E (I Aet{A{t))\-Ut)^oX'+{t) \ Z{t) = (^,0, 0)) . 
From Corollary [T] of Section |4] we know that 



E{\det{A{t))\.lA(t)^o 1 Z{t) = (u,0,0)) < u"" - 1 + 
Hence, 



'(0)f/^exp(-u^(24p"(0)-2)-^) 
v/24p"(0) - 2 



P{3T G S : X{T) = u, Thas minimal third coordinate on C„} 



(27r)3/2 



1 + 



''(0))=^''' exp (-^t^(24p"(0) - 2)-') 
x/24p"(0)-2 



(15) 
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• Secondly, if T is in the interior of a face Si , then, in this face, we choose the base ~^ such that 
is in the horizontal plane t\t2 such that along this vector, the second coordinate is not decreasing. 
Let us denote vector ~^ in the horizontal plane that is perpendicular to a and goes into S. It is easy 
to see that 

{X{T) = u, X(T) = 0, X'p{T) > 0, X\{T) < 0, X'^{T) < 0}. 
Apply Markov inequality and Rice formula to the field Y{t) = {X{t),X'a{t)), 

P{3 T & Si : X{T) = u, Thas the minimal third ordinate on Cu} 

< P{3teSi: X{t) = u, xUt) = 0, X'^{t) > 0, Xi,{t) < 0, X' < 0} 

< E(card{t e Si : X{t) = n, X'^{t) = 0, X'p{t) > 0, X'^{t) < 0, X^' < 0}) 
= E(|det(Y'(t))|Ix,(,)>oIx;w<oIx;'w<o I Y{t) = (^,0)) PY^t){n,0)dt 



'Si 

a2{ Si)(p{u) ^ 
/2tt 



|Xr(t)|X;,+ (t)Ix;w<o I Yit) = (u,0)) 



As in Theorem [T] it is clear that 



E(|XrWl I Y{t) = KO)) = v/l2p"(0) - 1^ 



E(x;,+ WIx;w<oiyW = KO)) 



Vl2p"(0) - 1 
1 - cos(/?, 7) 



+ 11$ 



vi2p"(o) - 1 y ' 



2V27r 

Observe that the angle between j3 and 7 is the angle 61 between the face Si and the horizontal plane, 
then the probability that there exists one point with minimal third coordinate on the level set and 
in the interior of the face Si is at most equal to 



0-2{Sl)(p{u){l — COS 6*1) 

47r 



\/l2p"(0) - lip 



+ u$ 



Vl2p"(0)-1 



Taking the sum of all the bounds at each faces, observing that 

n 

^ (72(50 COS = 0, 

i=l 

we have the following upper bound for the probability of having a point T with minimal third 
coordinate on the level set and belonging to a face: 



a2(S)(p{u) 
An 



\/l2p"(0) - Ip 



v/l2p"(0) 



(16) 



V12P"(0)- 

• Thirdly, when T belongs to one edge, for example Fi. Let us define 'rf is the upward direction 
on this edge, i.e such that along this vector, the third coordinate is not decreasing, and 7^ and ~^ 
are the two horizontal directions that go inside two faces containing the edge. Then, 

{X{T) = u, X;(T) > 0, XUT) < 0, X'^{T) < 0}. 

By Rice formula, the expectation of the number of the points in Fi satisfying this condition is 

J E (X+(t)Ixi(t)<oIx^(t)<o I X{t) =y}j X px(t)(u) dt 
= ai(Fi) ^(u) E (^X'+{t)Ix'^it)<olx'^(t)<o) ■ 

Let it and b be two vectors in two faces containing the edge Fi and perpendicular to if ; 61 be the 
angle between it and if ; 82 he the angle between ~^ and ~^ . It is clear that 



X (i) = cos 6I1 (f ) + sin 9i X (t) , 



12 



X'^it) = cos6»2 X'^{t) + sin 6*2 
and cov{X'a{t), X'i,{t)) = cos 63, where 63 is the angle between two faces containing the edge Fi. Then, 

E (^X'^'^lx'^(t)<o^x'^(t)<oJ 

=E ^X^+I^cos $1 (t)+sin $1 X'^ (t)<0}I{cos $2 X'^ (t)+sin 62 Xl(^t)<0}^ 

POO 

= X ifiix) F(x) dx, 
Jo 



where 



So, 



F{x) -E (ji-{co3ei x^(t)+sinei x'^(t)<o}i{cose2 x^(t)+siiie2 xl(t)<o} I X'nit) — 



. I ..^ /-cot^i.a; + cos03 cot6'2.a; 

F {x)= — cot Qi cot 62 -x) * : — 

\ sm^s 

^„ . _^ /- cot ^2 -a; + COS03 cot^i .a;\ 

- cot 9i (p{- cot(6'i x)) $ — , 

\ sm ^3 / 

By integration by parts, 

POO POO 

I X ifi{x) F{x) dx = - F{x) d{(p{x)) 
Jo Jo 

= F{0)ip(0)+ / ip(x)F'{x)dx. 
Jo 

It is easy to check that 

1'°° /SI/ s , 1 arctan(m) 
ifi{x)^{mx)dx = - + 

r / \ J 1 arctan(m) 
J ^{x) $(mx) dx= . 

Prom the above results, we have 

F(0)^(0) = ^ /" ^(2/)$ (^^^y) dy 



and 



f°° / \ TP'/ \ J f°° f \ I i-Q . ^ f -cot 61.x + 008 03 COt92-x\ , 

/ (p(x) P [x) dx = — cot 6*2 / ip(x)ip{— cot 6*2 -a;} "P -. — dx 

Jo Jo \ smBs J 

f°° r \ I x ^ / - COt 6^2 + COS 6*3 COt Q\ .x\ , 

— cot^i / (p(x)ip(— cot ffi .X) ^ \ : — 2i I dx 

Jo \ sin6'3 / 

— cos 02 / 1 1 / — sin $2 ■ cot 9i + cos 63 cos 62 
= — 7 + arctan — 

1/277 V4 27r V sin6l3 



— cos( 



'1/1 1 / — sin 61 cot 62 + cos O3 cos 9\ 

— T + 7^ arctan ^ — 

V4 27r V sin6l3 



Therefore, the probability that there exists one point with minimal third coordinate on the level set 
C{u) and belonging to F\ is at most equal to 

cri(-Fi)(7r — 63 )(/?(«) [—00362 /l 1 / — sin 02 cot 61 + cos 03 cos 62 
/o^^3/2 + '^i(^i) 7 + ^ a-^ctan — 



+ 



V27r V4 27r V sin6l3 y 

— cos 9\ f 1 1 / — sin 01 cot 02 + cos 03 cos 0i 

7 + arctan — 

^/2^ V4 27r V sm03 
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Summing up all the terms at all the edges, we obtain the bound 



(27r)3/2 



+ ip{u)Y,Gt[F,) 



COS 62i I I 1 / — sni 621 cot On + cos Osi cos 621 

H arctan 

sm tl-ii 



V4 27r 



— cos 6ii f I 1 / — sin 9u cot 62i + cos Osi cos ^1, 

+ —= — — I arctan 

V 4 27r V am 63, 



By definition, 
Now, we prove 

n 



^ai(F0(7r-6»30 =47rA(5'). 



cos6'2i fl 1 , / — sin 6*2^ cot 6'ii + cos Ssi 008 6*21 
— — 1 arctan 

sin dii cot 6*2^ + cos Osi cos Si, 



cos t'li / 1 1 
+ — T + arctan 



y/2n \4: 2tt \ smt/3i 

Indeed, from Lemma [l] we have 

— sin 9\ cot 62 + cos 63 cos ^1 — cos h 
sin S3 sin h ' 

where /i is the dihedral angle at "c^, i.e, the angle between the horizontal plane and the face containing 
1^ and 1/. Since h is constant for each face. 



/= ^ ? cos 6*1 f - + — arctan | ^ ) ) = 0. 

se{Si,...,Sfc} !cs ^ ^ 



sin /i 



Therefore, we have the following upper bound for the probability of having a point T with minimal 
third coordinate on the level set and belonging to an edge: 

2A(S)93(m) 



(27r)i/2 



From (151, (I61, (171 and the fact that 'P{X{0) > u} ^ ^{u), the result follows. 



(17) 

□ 



4 Computation of the absolute value of the determinant 
of the Hessian matrices 

As we see in the proof of Theorem |2j we deal with the following 

E(| Aet{X"{t))\lx..(t)^o I X{t) = u, X[{t) = 0, X^(t) = 0). 

To evaluate this quantity, we have the following statement that is one of our main results in this 
paper: 

Theorem 3. Let X he a standard stationary isotropic centered two-dimensional Gaussian field. One 
has 



E{\det{X"m I {X,X[,X^){t) = {u,0,0)) = u^-l + 2^- 



-"{0)f^^ cxp(-u^(24p"(0) - 2)-^) 



(18) 



v/24p"(0) - 2 

Proof. Under the condition, the vector (X"i, ^"2, ^22) has the same distribution with (11,^2,^3) + 
(—It, 0, — u), where {Yi, Y2, Y3) is a centered Gaussian vector with the covariance matrix: 



E = 



/ 12p"(0)-l 4p"(0)-l \ 

4p"(0) 

y 4p"(0) - 1 12p"(0) - 1 y 
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Then, the LHS in ( 18 1 can be written as 



E{\X^,{t)X^,{t) - XUtf\ 1 (X,X(,X0(t) = KO,O)) 
E{\{Y^~u){Y3-u)-Yi\) 
E{\YiYs - Yi - u{Yi + Ys) + u^\) 
E(\ <Y,AY > + <b,Y > 





( 





1 






( 




where A — 







-1 




and b = 









[ 


1 

2 





J 




[ 


-u ) 



Here, from Theorem 2.1 of [11) . the expectation is equal to 

2 



E(| < F, yiy > + < 6, y > +u^\) = 



r^{l~ F{t) - F{t))dt, 



where 



It is clear that 



and 



F(t) = 
F{t) 

m 



ex.p(itu^ - 2"4^ <b,{I - 2itT.A)-^T,b >) 
2det(J-2jtEyl)i/2 

exp{itu^[l - it{16p"{0) -2)]-'^) 
^ 2(1 + Sitp"{0))[l - if(16p"(0) - 2)]V2 ' 

exp(-itu2[l + it(16p"(0) - 2)]"^) 



2(1 - 8it/9"(0))[l + it(16p"(0) - 2)]i/2 



F(-f). 



So, the expectation is equal to 
1 



t' 



Al - Fit) - F(t))dt = Re(- 
= Re(- 



(1 - 2.F{i))dt) 

exp(itu^[l-it(16p"(0)-2)]"^ 



(1 + mp"{0))[l - it(16p"(0) - 2)]'/ 



Here, we apply the residue theorem to compute 



1 _ exp(it.."[l~»f(16/(0)-2)]-^) 



(1 + 8itp"(0))[l - it(16p"(0) - 2)]!/^ 
= 2i. (sum of residues in upper half plane) + i- (sum of residues on x-axis) 

The residues come from two poles at i.(8p"(0))^^ and and we see that: 
The residue at is equal to 

exp(ite2[l - jt(16p"(0) - 2)]-^) 



dt 



-i.u + i. 



t=i.{8p"(0))-l 



(1 + 8jfp"(0))[l - ii(16p"(0) - 2)]i/2 

And the residue at i.(8p"(0))^^ is equal to 

(1 + 8^^p"(0))[l - it(16p"(0) - 2)] - exp(aM^[l - a(16p"(0) - 2)]"^; 

t2.8ip"(0).[l - jt(16p"(0) - 2)]i/2 
^ (8p"(0))«/^exp(-^.^(24p"(0)-2)-^) 
\/24p"(0) - 2.i 

These two residues imply the result. □ 
We have the corollary 

Corollary 1. Let X be a standard stationary isotropic centered Gaussian field. One has 

(8p"(0))3/^ exp(-K^(24p"(0) - 2)-^) 



E(|det(X"(f))|.Ix"(t)^o I {X,X[,X'^)(t) = («,0,0)) <u2-l + ^ 
Proof. The result follows from two observations 



v/24p"(0) - 2 



15 



. |dct(X"(t))|.Ix"(t)^o < ^ 



det{X"{t))\+det{X"{t)) 



. E(det(X"(t)) I {X,Xi,X'2m = (m,0,0)) = - 1. 

□ 

5 Numerical comparison 

In this section, we compare the upper bounds given by the direct method and record method 
with the approximation given by the EC method. For simphcity we limit our attention to the case 
where S is the square [0,T]2 and X is a standard stationary isotropic centered Gaussian field with 
covariance function p{\\s — Note that only p"(0) plays a role, the exact form of p does not need 
to be specified. More precisely, we consider 

1. the approximation given by the EC method 

Pe{u) = #(m) + ^==V5(m) + 75- 
V27r ^TT 



2. and the upper bound given by the direct method 

Pm{u) = *(w) + / [cip{x/c) + x^{x/c)] (p{x)dx 

v27r Ju 



2tt 



<p{x/c) 



ip{x)dx, 



where c = -/l2p"(0) - 1, 



3. and the one given by the record method 

¥r{u) = $(m) + -^=ip{u) + — [ap{u/c) + M*(w/c)] </?(«). 
V27r 27r 

It is easy to see that Pb is always less than Ph and Pm- We will prove that Pr{u) is smaller than 
Pm(w) as u is large. Indeed, if we compare the "dimension 1 terms" (corresponding to cti(9S')), we 
have 

[c^p{x/c) + x^{x/c)\ ifi{x)dx — ip{u) 
= !u [cv(a;/c) + x^{x/c)\ (p{x)dx - xif{x)dx 
= !u [c¥'(a;/c) - x^{x/c)\ ip{x)dx > 0, 

since when x > 0, 

^>^ix). 

X 

So the term in the direct method is always larger when m > 0. 

Let us consider now the two terms corresponding to 1X2 (S): 

3/2 



• Ar = [ap(u/c) + w$(w/c)] ip{u) = U(p{u) + A 
It is easy to show that , as m — ^ +00, 



ifi{x)dx = uip{u) + Ad. 



j Lp(^^^ip{x)dx = {const)<b^u ^^'^'^ 



and that 



Ar ~ {cons 



{const)u ^'p\u 



^2 
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This shows that for u sufficiently large Ar is smaller than Ad. 

The numerical comparison is performed in Figure 4 for six different situations, ft shows that the 
record method is always better than the direct method. EC method and record method are very 
close, but it is not possible to identify the better among those two since Pe can be smaller than the 
true value. 



1.2 




"0 5 1 1 5 



□<u<1 5 

(a) T = 1, p"{0) = 0.25 




□<u<1 5 

(b) T = 1, p"(0) = 0.5 
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